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Abstract 

Survival in groups of termites is known to be socially facilitated 
since individual longevity is increased together with the size of the 
group. Here we report on experimental evidence that survival as a 
function of group size in Cornitermes cumulans (Kollar) does not in- 
crease asysmptotically but it peaks in a group size range at which 
individual longevity is maximal. Strinkly, this same effect in other 
social insects was noted more than fifty years ago, but it was regarded 
as "aberrant" since it has been generally assumed that individual sur- 
vival should increase with the size of the group. We also provide 
evidence that individual activity -measured as mobility- is maximal 
at approximately the same group-size range as the observed survival. 
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1 Introduction 



Grasse and ChauvinI (|194J) reported a number of experiments to check the 
effects of group size on the survival of Apis sp., Polistes sp, Leptotorax 
sp., Formica sp. (Insecta: Hymenoptera) and Reticulitermes sp.(Insecta: 
Isoptera). In all of them, but in Polistes, survival increased with increasing 
group size. In Polistes they reported an "aberrant" result in which survival 
in both large and small groups was smaller than in the middle ones. The 
authors atributed these obeservations to bad experimental procedures or to 
a very particular physiological trait of the species. The experiment was hence 
repeated several times -even by independent observers- but the "aberrant" 
result was always found. On a survey experiment performed on the sum- 
mer 1994, we test grouped individuals of Syntermes sp., Nasutitermes sp., 
and Cornitermes sp. (Insecta, Isoptera). Preliminary results showed the 
same unexpected behaviour on the survival of grouped Cornitermes cumu- 
lans (Kollar). 

Such a convergence of results allow us to suspect that Grasse & Chauvin 
(1944) have in fact unwarrantly spotted a pattern that is biologically sound. 
We base our suspicion on previous experimental works which show that (i) 
individual survivorship of termites, while affecte d bv group size, mav be re- 



lated to individual mobility of group components (jMiramontes and DeSouza , 




19961 ). and (ii) longevity of individual termites seems to obey a hump-shaped 
curve, clearly depend ent on some (yet) unknown attribute of group size 



(DeSouza et al 



20011 ). In addition, there is theoretical evidence that co- 



operation could either stabilize or destabilize the dynamics of a social group. 



3 



I 

leadin g to maximal individual fitness at intermediate colony sizes (jAviles, 

This paper aims to present experimental evidence linking individual mo- 
bility, survivorship, and group size in termites. To achieve this we show that 
both individual mobility and survival are humped functions peaking at sim- 
iliar group-size range. We argue, on experimental and theoretical grounds, 
that survival of starving grouped termite workers -at least in Cornitermes 
sp- obeys a humped-shape curve which present a biologically sound, not 
fortuitous, relationship to individual mobility. 

2 Material & Methods 

To confirm the trend observed in survey experiments, two further experi- 
ments were performed in 1996 and 1997 with grouped termites of this species 
collected from wild colonies in Vigosa in the state of Minas Gerais, South- 
eastern Brazil. The first experiment, spanning several days, was designed to 
measure longevity while the second was intended to measure individual ac- 
tivity. Both experiments involved termite workers (third instar and beyond) 
of different group sizes. 

In the longevity experiment, termites from a wild colony were randomly 
placed in groups of 1, 2, 4, 8, 12, 16, 20, and 24 (each group replicated eight 
times). Collection took place in 04 April 1996. The groups were confined in 
test tubes made of transparent glass (9.5 x 1.4cm) with hermetically sealing 
rubber caps. Each test tube is a 'replicate'. Test tubes have been previously 
washed and rinsed, soaked in sodium hypochlorin for 24h, rinsed again, and 
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sterilised at 180 °C for one hour. Tubes with termites were kept horizon- 
tally, separated by plastic foam to prevent stridulation or other mechanically 
transmitted signals to propagate between the tubes. The workers were incu- 
bated in the dark at a constant temperature chamber (25 °C ± 0.5) and were 
allowed to acclimatise for 12h. No food or water was provided, but tubes 
were opened each 24h, to allow air exchange. Tubes were exposed to light 
during the counting of survivors only (no more than 5 min). Observations 
were made each eight hours intervals, until all individuals were dead (153h; 
20 observations). Termite groups in which dead individuals presented any 
sign of cannibalism were excluded from data analysis. 

In the activity experiment, termites from another wild colony were ran- 
domly placed in groups of 2, 4, 6, 8, 16, 20, and 30 individuals (each group 
replicated four times). Collection took place in 12 October 1997. The exper- 
iment was setup as above, using a controlled temperature room (24 °C ± 0.5) 
rather than a chamber. Termites were allowed to acclimatise in such a room, 
for three hours before the observations began. Two observers took notes 
of the number of termites moving in each group, at five minutes intervals, 
during ten continuous hours. Each observer managed two replicates of the 
experiment (two tubes of each group size). Light intensity in the room was 
the lowest possible to allow the work of the observers. 

Statistical analyses aimed to check whether group size would affect (i) 
survival and (ii) mobility of termites. Either quantities were treated as the 
y-variable in separate regression analyses, always with group size as the x- 
variable. Survival was measured as the average number of hours spent until 
individual termites die in a given group, which was calculated for each repli- 
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cate (a given tes t tube containing termi t es'), u sing Weibull frequency distri- 



bution (^Crawley 



Pinder III et al, 



19781 ). Survival functions obey the 



general form: 



In S{t) = li-^'f 



(1) 



where S{t) is the proportion of individuals from the initial cohort (the initial 
group of termites) still alive at time t, a is the shape parameter, /i is the 
mean number of hours until death by the termites from a given replicate. 
The first step in the survival analysis was to establish the value of the shape 
parameter a, thereby defining whether or not the proportion of individuals 
dying was constant through time (a = 1, dN/dt = l/fi, exponential fit). To 
do so, survival curves were tentatively fit to the data, simulating different 
shape parameters, one at a time, until the error deviance reached an asymp- 
tote. The model obtained was then compared with a model assuming a = 1 
(exponential fit), on the basis of their contribution to the error deviance. 
Parsimony requires to choose the exponential fit (simplest model) if it does 
not differ significantly (P > 0.05; tables; 1 df) from an alternative, more 
complete, model. 

After defining the shape parameter, the next step was to calculate the 
mean number of hours termites spent to die when confined in each group size 
(/i in equation 1 above). Each of such means represent the number of hours 
spent to die averaged across all individuals belonging to a given replicate, 
and are referred to as "mean time to death". Mean times to death for each 
replicate were then collapsed into a single arithimetic mea n for each grou p 



size, a valid procedure to avoid pseudoreplication effects (jCrawlevl . 
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Such a procedure produces eight means (each one represented by a dot in 
Fig.l), one for each group size. 

Means thereby produced were used to check whether termite survival (y- 
var) would be affected by group size (x-var). To do so, models were fitted 
to the data, starting from a null model {y = b, were b is the grand mean), 
and adding new terms until achieving the best trade-off between percent of 
variance explained (r^), and P- values obtained (see Table 

The general pattern of mobility of termites in a given replicate (a test 
tube containing termites) was measured by fitting a simple regression line of 
the form y = b + ax, through the data points formed by plotting the number 
of moving termites (y-var) against time in minutes (x-var). The slope a 
of each line was used as a measurement of the general pattern of mobility 
for termites belonging to that replicate. If the number of moving termites 
increases as times goes by, the curve fitted to the data would show a positive 
slope. Conversely, a negative slope would indicate a decay in group mobility. 
When no clear trend is to be observed for a given replicate, its "mobility 
curve" would show a null slope. For each group size we obtained four slopes, 
since there were four replicates (i.e., four test tubes containing termites). The 
slope values thereby obtained were collapsed into a single arithimetic mean 
for each group size aiming, as above, to avoid pseudoreplication effects. Such 
a procedure produces seven means (each one represented by a dot in Fig.2), 
one for each group size. Such means are refereed to as the "general pattern 
of mobility" for that group size. 

Means thereby produced were used to check whether the general pattern 
of mobility (slope a of mobility lines, see above) would depend on the size of 
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the group were termites have been confined. Models were fitted and selected 
as above (see Table |2)). 



All statistical analyses were performed using R (jlhaka and Gentleman 



19961 ). 



3 Results 

The mean time to death of starved termites obeys a humped function of 
group size, presenting a peak at a characteristic group size and decreasing 
at both smaller and larger densities (Polynomial model, order four: F[4;3] 
= 20.22 ; r'^adj = 0.916; P = 0.02; Fig. 1; Table ^ Accordingly, group 
activity presents a similar dependency on group size, showing a peak at 
an intermediate density and decreasing at both smaller and larger values 
(Quadratic model: F[2;4] = 9.095; r^adj = 0.730; P = 0.03; Fig. 2; Table 
121). This quadratic model explains the data better than the simple linear 
one: adding a quadratic term to the linear model, increases by 18.7% the 
percentage of variance explained and improves the P- value (Table |2)). 

Such results are consistent across different procedures of data analysis. 
As specified above, each dot in Figs.l and 2 represents an arithmetic mean 
taken across the respective replicates. If, however, we plot the raw data (i.e., 
splitting the arithmetic means into its components), we observe the same 
general shape in both curves (Figs. 3 and 4). 
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4 Discussion 



Across the Isoptera, a range of different socially facilitated behaviours has 
been extensively reported, although survival is the one that has attracted 
more attention (Table El). In all survival studies, longevity has been reported 
to increase as group size increases. This may lead one to intuitively assume 
an asymptotic behaviour of such a curve, accepting that very large groups 
would attain the maximal possible survival rates. In this sense, the results 
reported by Grasse & Chauvin(1944) in which survival of wasps is maximal 
at moderate densities (rather than at large group sizes), could be consid- 
ered "aberrant". However, a striking similarity can be observed between 
such results and those found here for C. cumulans termites (Fig. 1). More 
puzziling, termites from different wild colonies of the same species show a 



similar humpec 



( DeSouza et al 



'unct ion, even when they have been poisoned by insecticide 



20011). Such a convergence of results lead us to suspect that 
this pattern is, in fact, biologically sound. It seems that survival of starv- 
ing termite workers -at least in C. cumulans- is not asymptotically related 
to group size, but obeys a humped function, in which an optimal density 
(=number of individuals per unit area) assures maximal survival. 

Optimal densities in social insects have already been predict ed as a conse- 



quen ce of evolutionary pathways maximising individual fittness fjHigashi and Yamamura 



19931 ) . Accordingly, extraordinary lifespan is argued to d erive from evolution- 



ary p rocesses in which such insects benefit from sociality (|Keller and Genouc 



19771 ). Proximate mechanisms through which optimal densities would lead 



to better survival, however, are still poorly understood. Many works deal- 
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ing with cluster size in group organisms are based on optimality arguments 
regarding resou rce usage, territory defence, anti-predatory behaviour and dis- 



ease resistance 



1999; 



Robert! 



Creel 



199 



1997 



Fritz and Garine-Wichatitskv 



Rosengaus et al 



1996 



Giraldeau and Beauchamp 



1998^ . Others, inspired on novel con- 



cepts of the sciences of complex systems, would invo ke a relationship between 



group-size, infor r aation fluxes and task per 



Bonabeau et al 



Deneubourg et al. 



brmance (lAd 



mm 



Cordon. 



er and GordonL 



1992; 



Pacala et al 



19961 ). The interconnected nature of such systems is thought to be the basis 



for the self-organisation of a variety of cluster-re 



JM3 



Camazine 



1991 



Miramontes et al, 



ated phenomen ( 


Bonabeau et al. 


Theraulaz et al.. 


1995 


). Specif- 



ically f or termites, evidences seera to po int out that survival is related to mo- 



bility (jMiramontes and DeSouza . 



19961 ). In addtion, studies addressing in- 



terindividual interactions and self-orgnisation predict the existence of group- 
size effects causing info rmation transfer and behavioural diversity to reach 



near optimal conditions (jMiramonte 



Sole and Miramontes 



mM). Ac- 



cordingly, termites studied here present maximal values for group activity at 
intermediate densities (Fig. 2), which lie in the same group-size range as the 
observed survival (Fig.l). It seems,therefore, that the association between 
group-size, non-asymptotic survival and greater individual activity is not 
fortuitous in C.cumulans. Moreover, to the extent that mobility affects the 
rate of soci al contacts -a nd thereby, cooperation- our results seem to agree 



to those of lAvilesI (|l999t ). who show an explicit link between cooperation, 
individual fitness, and group size in social organisms. 

These results point out to specific mechanisms present in the dynamics of 
groups that may act to regulate observables such as density, number of indi- 
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viduals and rate of social contacts, that in turn have an impact on collective 
task performance. Arguably, behaviours such as colony fission, swarming, 
cannibalism, or resting, may be the natural mechanisms responsible for the 
regulation of the proportion of individuals in the colony. In other words, 
the patterns presented here may have unexpected consequences for under- 
standing the life cycles of these organisms, and therefore are worth exploring 
further. 
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6 Legend to figures 

Figure 1 Average number of hours termite workers spent to die when con- 
fined with conspecifics in test tubes, in the absence of food. Each dot 
represents an arithmetic mean across eight rephcates (eight different 
test tubes). Average number of hours for each replicate was calculated 
by censored survival analysis, with WeibuU distribution. 

Figure 2 General pattern of mobility of termites confined with conspecifics 
in test tubes, in the absence of food. Mobihty patterns are characterized 
here using the slope of a line fitted across a scatterplot of the number 
of moving termites (y- var) versus time (x-var) , four each of the group 
sizes tested. Positive slopes denote that the number of moving termites 
increases with time, negative slopes denote the opposite. Null slopes 
denote a constant number of moving termites across the time range 
when observations were done. Each dot represents the average of the 
slopes presented by each of four rephcates (four different test tubes) . 

Figure 3 Average number of hours termite workers spent to die when con- 
fined with conspecifics in test tubes, in the absence of food. Each 
dot represents the average number of hours termites survived in a sin- 
gle replicate (a test tube), as calculated by censored survival analysis, 
with WeibuU distribution. Within each group size, datapoints of this 
figure were collapsed into a single arithmetic mean in order to produce 
Fig.l. 

Figure 4 General pattern of mobility of termites confined with conspecifics 
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in test tubes, in the absence of food. Mobility patterns are characterized 
here using the slope of a line fitted across a scatterplot of the number 
of moving termites (y-var) versus time (x-var), four each of the group 
sizes tested. Each dot represents the slope presented by a single of four 
replicates (four different test tubes). Within each group size, datapoints 
of this figure were collapsed into a single arithmetic mean in order to 
produce Fig. 2. 
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mobility (slope of activity curve) 



-0.030 -0.025 -0.020 -0.015 -0.010 -0.005 0.000 

I I I I I I I 
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Table 1: Models explaining survival the effect of group size (x-var) on the 
mean time to death (y-var) of starved termites confined in test tubes. Models 
provided here are intended to be illustrative of the general shape of the curves 
described by the data, rather than a statement on its specific poisition when 
plotted. Models were calculated on the datapoints of Fig.l. Polynomial 
model of order 5 is not significant. 



Model 


r^adj 


P 


y = 66.829 




1.00 


y = 88.8355 - 15638X 


0.309 


0.09 


y = 75.8990 + 0.8722a; - O.lOie^^ 


0.262 


0.20 


y = 47.94922 + 16.77077x - 1.74097^2 + 0.04376x3 


0.798 


0.02 


y = 29.783398 + 31.566777x - 4.431666x2 + 0.211764x=^ - 0.003355x^ 


0.916 


0.02 
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Table 2: Models explaining survival the effect of group size (x-var) on the 
mobility (y-var) of starved termites confined in test tubes. Models provided 
here are intended to be illustrative of the general shape of the curves described 
by the data, rather than a statement on its specific poisition when plotted. 
Models were calculated on datapoints of Fig.2. Polynomial model of order 3 
is not significant. 



Model 



r^adj P 



y = -0.009973 



1.00 



y = 0.0008778 - 0.0008832a; 



0.543 0.04 



y = -8.299 -10-3 + 1. 110 -lO-^o;- 6.396- 10" V 0.730 0.03 
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Table 3: Some different behavioural traits that are known to be socially 
facilitated across the Isoptera 





± 1 dl u 




JDctilCOoliof fflCh 


V^VdllUIIl Qt;Vt:jiOpiIlt;IlIj 


Vjidsse 1 iyoy i 


As Tl H IT T'TYl O Q 

JD LJ LUjliC.1 iiLt'O 






LytpilUlUbtl IlLto 


V^VcLIlUlll U.t;Vt:lUpillt;llL 


vjidbse 1 ±yoy 1 


CoptoterTnes 


Survival and feeding 


L/Cnz ana vvunanis iiyoui 




loieiduce Lo puisoning 




Cvyptoterfnes 


Survival 


vv imams ex ai. i lyozi 


Kalotermes 


Nest digging 


Snringhetti (1990) 


Macrotermes 


Caste differentiation 


Okot-Kotber n988l 


Nasutitermes 


Survival and feeding 


Lenz and WiUiams flQSO) 




Survival under starvation 


Miramontes and DeSouza flQQG") 


Reticulitermes 


Survival 


Grasse and Chauvin fl944~) 


Zootermopsis 


Survival under infection 


Roseneaus et al. fl998) 
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